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Comparison between theoretical four-loop predictions and Monte Carlo 
calculations in the two-dimensional N- vector model for N = 3, 4, 8 
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We have computed the four-loop contribution to the beta-function and to the anomalous dimension of the 
field for the two-dimensional lattice TV-vector model. This allows the determination of the second perturbative 
correction to various long-distance quantities like the correlation lengths and the susceptibilities. We compare 
these predictions with new Monte Carlo data for N = 3, 4, 8. From these data we also extract the values of various 
universal nonperturbative constants, which we compare with the predictions of the 1/N expansion. 



The A- vector model describes configurations of 
classical spins taking values on the unit sphere 
gN-i [fjA^ -yVe consider here the standard 
nearest-neighbor action 

n(o-) = -f3j2<T.-<Ty ■ (1) 

The perturbative renormalization group predicts 
that when [3 / oo, the (infinite- volume) long- 
distance quantities of lattice theory should give 
the same results as the continuum theory in the 
MS normalization, provided that one rescales 
lengths by the factor 
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Here the constants C# are universal but non- 
perturbative, while the coefficients a;, 6; and di 
can be determined at the (i + 2)-th order of the 
perturbative expansion. There is a prediction [1] 
for the exponential correlation length (= inverse 
mass gap) in the isovector channel: 



For the isovector and isotensor two-point func- 
tions 
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As the model should not have bound states, in 
the isotensor channel one expects 



we shall consider the RG predictions for the cor- 
relation lengths 
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The constants related to the second-moment 
correlation lengths and the susceptibilities are 
known analytically only in the large- A expansion, 
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namely, through order l/N both in the isovec- 
tor [2] and isotensor [3] sectors. In particular, 
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so that even for N = 3 this ratio differs only 
marginally from 1 (in good agreement with Monte 
Carlo simulations [4]). 

We can also try an improved expansion param- 
eter [5,6] based on the isovector energy Ey = 
(fo,o • 0"i,o): we define 
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This f3 e g has the property that at N / oo with 
[3 = [3/ N fixed, there are only exponentially small 
corrections to the two-loop predictions for corre- 
lation lengths and susceptibilities. 

A precise confirmation from Monte Carlo sim- 
ulations of (5) with the correct nonperturbative 
constant (8) / (10) would be good evidence in fa- 
vor of the conventional asymptotic-freedom pic- 
ture [7], which has been criticized [8]. At the 
last Lattice conference [9] we presented data for 
the case N = 3 at infinite-volume correlation 
lengths up to k, 10 5 [10], obtained by using 
finite-size-scaling extrapolation at fixed [3 [11]. 
The discrepancy of these data from the three- 
loop predictions [12,6] was already quite small 
(~ 4%), and in the "improved expansion parame- 
ter" the agreement was even better. For this rea- 
son we considered it worthwhile to compute the 
next-order perturbative correction, in order to see 
whether the remaining discrepancy — which is 
nevertheless larger than the estimated statistical 
error — can be removed. We therefore computed 
the four-loop contributions to the beta-function 
and the anomalous dimension of the field for the 
lattice model (1); from these we derived the 1//3 2 
corrections to the correlation lengths and the gen- 
eral spin-n susceptibility [13]. For example, we 
have 
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In Figure 1 we plot ^ estimate divided by 
the two-loop, three-loop and four-loop predictions 



f i l I I 

i 
I 



i t 



2.0 



2.5 



3.0 



Figure 1. ^estimate! ^[theor VCTSUS I 3 for tlle 

0(3) model. Error bars are one standard devi- 
ation (statistical error only). There are five ver- 

(2) 

sions of ^ theor - standard perturbation theory in 
1//3 gives points + (2-loop), x (3-loop) and hJh (4- 
loop); "improved" perturbation theory in l//3 ej ff 
gives points □ (2-loop) and O (3-loop). 
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(points +, x and hJh) or by the "improved" two- 
loop and three-loop predictions (points □ and O). 
The four-loop truncation of (5) is now fully com- 
patible with our last extrapolated point. It falls 
roughly halfway between the three-loop trunca- 
tions in 1//3 and 1 / f3 e g. 

A similar result is obtained for N = 4, where 
the central estimate of our last point (which is at 
the rather small correlation length ^ « 155) dif- 
fers from the theoretical four-loop prediction by 
only 4-5%; and for N = 8 [14], where the central 
estimate of our last point (which is at ^ « 650) 
shows an extremely good agreement (better than 
1%). 

We have also tried to take into account higher- 
loop corrections by using information from the 
1/A-expansion. Let 

a^l N \N) = N n ~ 1 a n (13) 

be the leading contribution to a n in the limit N / 
co. We computed a n up to n = 8 [13]. How good 
is the approximation in which only such a term 
is retained? We can compare with the known 
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coefficients a\ and 02. For N = 4, we have 
«i( 4 ) _ o -70. «2( 4 ) 



3.73; 



«)' J (4) a 2 1/JV; (4) 
while for N = 8, we have 
«i(8) 1 m «2(8) 



1.91; 



Ul '(8) 4^(8) 

The convergence seems slow; indeed, only at N > 
50 are the first two coefficients correct to within 
10%. Let us now define k n by 

a„(8) = *„8 n_1 an . (16) 

Already at /? = 5.80 « 33) the Monte Carlo 
value [14,15] is in good agreement with the theo- 
retical predictions. Indeed we get 
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0.998 when k n = 1 Vn 
1.001 when k n =2Mn 
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with a statistical error of ±0.002. For larger val- 
ues of [3 this ratio remains roughly constant, al- 
though the error bars grow. We can thus claim a 
nice control of (5) for N = 8. 

Having verified (5), we can now extract from 
the Monte Carlo data a numerical evaluation of 
the nonperturbative universal constants C#. For 
the limiting ratio £y ^ we have the numeri- 
cal results 3.51 for N = 3, 3.14 for N = 4, and 
2.77 for N = 8, with error bars less than ±0.01, 
which can be compared with the l/7V-expansion 
prediction [3] 
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for N = 3 

for N = 4 (19) 
for N = 8 

To extrapolate the asymptotic value of the 
isovector susceptibility, we found it convenient to 
study the dimensionless ratio 
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because our knowledge of the lattice beta and 
gamma functions at four loops allows the determi- 
nation of the first three coefficients c; [6,13]. We 
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Figure 2. Estimate of C* Xv [from (20)/(8)/(10)] 
versus [3 for the 0(8) model. Error bars (one 
standard deviation, statistical error only) shown 
for clarity only on one set of points. There are 
seven versions of xv /(Cv ^) 2; standard perturba- 
tion theory in 1//3 gives points ± (leading), x 
(with ci), ± (ci )2 ) and ^ (01,2,3); "improved" 
perturbation theory in 1/ (3 e g gives points □ (lead- 
ing), O (ci) andjqj (ci 2 ). 



then use the exact formula (8) / (10) to estimate 
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see Figure 2 for N = 8. We find 



r 10.8 ±0.8 for # = 3 
C xv = I 5.9 ± 0.1 for N = 4 
I 5.6 ±0.1 for # = 8 



(21) 



which can be compared with the 1 /TV-expression 
[2] 
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3.67 for N 
4.32 for # 
5.30 for N 



where je is Euler's constant and 
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Clearly the 0(1/ N 2 ) corrections are significant! 
We have made the wild guess that the exact ex- 
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Figure 3. C Xv versus 1/N . The straight line is 
the 1 /TV-prediction (22), while the curve is the 
wild guess (25); the points O are our Monte Carlo 
estimates. 



pression for C Xv is 
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in analogy with (8). In Figure 3 we compare our 
Monte Carlo results with the 1 /TV-prediction (22) 
and the wild guess (25). The result for N = 3 
is close to the guessed formula but outside the 
statistical errors, while the values for # = 4,8 
follow the Ansatz nicely. It would be interesting 
to test (25) by computing the 0(1/ N 2 ) term in 
(22). ^ \ 

The agreement with the l/7V-expansion is 
much poorer for the isotensor susceptibility: our 
Monte Carlo data yield 

f 1200 ± 100 for N = 3 



C XT = < 23 ± 2 for N = 4 

I 4.7 ±0.2 for # = 8 
compared to the l/7V-expansion result [3] 
0.0296 
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Clearly the 0(1/ N 2 ) and higher corrections must 
have a drastic effect for N < 20. 
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